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Design of thick-walled components for power plants 
from crack-growth models 



DESCRIPTION 



Technical Field 

The invention relates to the field of power plant design. It 
relates in particular to an iterative method and a computer 
program product for designing components for industrial plants, 
in particular thick- walled components for power plants, as 
claimed in the precharacterizing clause of claim 1. 



Prior Art 

Thick- walled components appear frequently in power plants, for 
example as tubing for superheaters, reheaters, economizers and 
as coolers, drums, rotors, turbine and compressor blades and 
casings. Typical component design approaches take conditions 
such as static water/steam flows, static operating temperatures 
and static stresses into account in a steel selection. In one 
common design approach, thick-walled components for power 
plants are normally designed to last a given number of hours at 
a given full load. Full-load equivalent -operating-hours (FLEOH) 
for each component, e.g. in a boiler, are determined by means 
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of look-up tables which are widely available to persons skilled 
in the art. To find a component's FLEOH, both steady-state 
full-load temperature and steady-state full-load stress the 
component is subject to under steady-state full -load conditions 
are required^ as are material properties, e,g. steel type, of 
the component. When applying this approach, stronger materials 
will typically be chosen for tubes and welds in superheated 
parts of a boiler when designing the power plant, thus provid- 
ing such components with FLEOHs similar to tubes and welds in 
low temperature regions. The FLEOH number for a component cor- 
responds to a life expectancy of that component when constantly 
subjected to the steady-state full-load temperature and the 
steady-state full-load stress. If values for temperature and 
stress during operation exceed steady-state full-load values 
assumed when designing the component, the life expectancy of 
the component will decrease, or, in other words, ageing of the 
component will be accelerated. As a consequence, the component 
will fail earlier than expected. 

The FLEOH approach is conservative. In general, a power plant 
will not be operated at full-load at every given moment of its 
lifetime. An alternative approach is to use load- dependent 
equivalent-operating-hours (EOH) . The EOH is an average esti- 
mate of a component's life expectancy based on typical load 
profiles and typical numbers of start-up and shutdown sequences 
that occur during a components lifetime. As a result, a more 
accurate estimate of the life expectancy is achieved and a less 
conservative design is possible. 

In both approaches, however, to ensure that components do not 
fail unduly, harsher conditions than the ones that the compo- 
nents will finally be subjected to are assumed when designing 
said components. This leads to more or less over-dimensioned 
components, with, in general, thicker walls, larger diameters, 
materials and welds of higher quality, etc., than would actu- 
ally be required. This, in turn, leads to higher costs in 
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building corresponding plants, namely because of increased ma- 
terial consumption, higher transportation costs due to in- 
creased weight of the components, more difficult handling, etc. 



Description of the Invention 

It. is an objective of the invention presented herein to provide 
a method and a computer program product of the type mentioned 
initially that permits a design of components for industrial 
plants, in particular thick walled components for power plants, 
that avoids over-dimensioning. 

This object is achieved by the method and the computer program 
product as claimed in the independent claims. In an iteration 
method for designing a component, in particular a thick-walled 
component, of a power plant according to claim 1, a dynamic 
simulation of process variables is used as the basis for a 
crack-growth simulation. The process variables are computed by 
means of a process simulator for a given, time -dependent load 
profile. According to the invention, the load profile is time- 
dependent and the process simulator is capable of taking into 
account transient behaviour as it occurs e.g. during start-ups, 
shutdowns or load changes, so that the computed process vari- 
ables are also t ime- dependent . This has the advantage that com- 
ponents can be designed in such a way that their life expec- 
tancy under realistic operating conditions accurately matches a 
desired value. 

Further advantageous realizations can be found in the dependent 
claims . 



K6, 04.12.2002 



02/124 



- 4 - 

Brle£ Explanation o£ the Figures 

The invention will be explained in more detail in the following 
text with reference to exemplary realizations and in conjunc- 
tion with the figures, in which: 

Fig* 1 shows a potential location and geometry of ' a crack in a 
turbine blade. 

Fig. 2 shows a cross section of a pipe, with an example of a 
crack growth geometry • 

The reference signs used in the figures are explained in the 
list of reference signs. 



Detailed description of preferred eiobodixaents 

Starting point for the method according to the invention is an 
initial structure of a component, in particular a thick-walled 
component, said component forming part of a power plant, which 
in turn is characterized by an initial structure, wherein the 
initial structure of the component is comprised in the initial 
structure of the power plant. The word ^*structure" in this con- 
text refers to a physical, geometrical and material constitu- 
tion of the component and the power plant, respectively, and 
comprises shapes, dimensions, steel composition, weld materi- 
als, etc. For example, a combined cycle power plant comprises a 
gas turbine, a heat recovejcy steam generator, and a steam tur- 
bine. The steam turbine comprises a plurality of thick-walled 
components such as a rotor, a number of pipes, a casing, a num- 
ber of turbine blades, and a number of vanes. Each component is 
described by process dynamics, geometrical dimensions and mate- 
rial characteristics. In addition to the gas and steam turbines 
and the heat exchanger, a power plant may comprise further com- 



K6, 04.12.2002 



02/124 



- 5 - 

ponents, e.g. a compressor, a condenser, a feed water tank, a 
deaerator, coolers and several pumps. 

In what follows, the component need not necessarily consist of 
a single part, but may be a composite component comprising a 
number of parts that are joined together, e.g. by a weld. In 
this manner, joints like welds can easily be included in the 
design process. This is important since many failures encoun- 
tered in heavy-section piping occur at welded joints. Hence, 
not only the components base material, but also the weld re- 
gions require sufficient attention when designing power plants. 

In a first step, a response of the steam turbine and/ or the gas 
turbine to a t ime - dependent load profile is modelled by means 
of a process simulator. For the steam turbine, for example, 
elements such as rotor dynamics, steam temperatures, and pres- 
sure distributions within the steam turbine are taken into ac- 
count in a simulation. A corresponding modelling methodology is 
well-known to a person skilled in the art of power plant model- 
ling, and is, for example, described in detail in the book 
**Thermische Turbomaschinen" by Walter Traupel, published by 
Springer, Berlin, 1988, which will be referred to as Traupel in 
what follows. Implementation of the modelling methodology into 
a process simulator in turn is straightforward to a person 
skilled in the art of computer programming. Preferably, rotor 
dynamics, thermodynamic behaviour of pipes 4, and thermodynamic 
behaviour of steam within the turbine are taken into account in 
the simulation. From the process simulation, a plurality of 
process variables will result, comprising, but not limited to, 
densities, temperature distributions, pressures and/or veloci- 
ties of steam and/or fluids. For the gas turbine, pressure, 
density, velocity of the gas in a high and/or a low-pressure 
section of the gas turbine, pressure, density, enthalpy may, 
for example, be obtained by means of the process simulation. 
Temperatures of steam, water and air in a heat recovery steam 
generator may also be simulated. 
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The following mass balance equations are typically used to 
model the mass flow in a pipe. 

Vol = JJowr' - flow'Jr > 

where p is the density, p the pressure, and the remaining vari- 
ables are component specific constants. Similar equations can 
be written for energy balances and used for turbine and heat 
exchanger models. 

In a second step, stress exerted onto the component will be 
computed. The following types of stress are taken into account: 

• Tension stress due to centrifugal forces: pure tension due 
to the combination of rotor rotation and the blades own 
weight, for example computed using essentially the follow- 
ing formula: cr^(r*)= ^ f area(r)rdr, where a is the ten- 

area(r ) 

sion stress, CD a rotational velocity, and p a component 
density. 

• Bending stress due to pressure gradients, for example com- 
puted using a formula taking into account axial and longi- 
tudinal bending moments as described in Traupel. 

• Thermal stress due to temperature gradients within the 
blade, for example computed using essentially the follow- 

mg formula: 'ztT'^ — ^ + — r = ^/^(Tr-^x)' where ^ is the 

dr r dr dz ^-^^r 
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Stress potential, Tr the temperature, Tr* a mean tempera- 
ture, and z and r polar coordinates of the component. 

• Dynamic stress due to fast oscillations of process vari- 
ables caused by flow fluctuations and/or excited eigenfre- 
guencies. This type of stress is especially important dur- 
ing start-ups, shutdowns and load changes. For example 
pipe stresses are directly proportional to internal pres- 
sures that vary dynamically. Similarly as can be seen from 
the above mentioned stress formulas, thermal and bending 
stresses are directly related to process variables such as 
temperature and pressure that also vary dynamically. 

Computational procedures for modeling the various stress types 
are again known to a person skilled in the art of power plant 
modelling, and are once more described in Traupel. 

In a third step, a life expectancy of the component is deter- 
mined for the given load profile. For this purpose, a hypo- 
thetical crack is introduced in the component, and a growth 
with time of said hypothetical crack is modelled by means of a 
crack-growth simulator. The hypothetical crack reflects the 
fact that, even when out -of -the -box new, the component has a 
number of microscopic or minimal cracks. In addition, it is 
further assumed that each component has a critical crack 
length. The life expectancy is the time it takes for the hypo- 
thetical crack to reach the critical length. Depending on the 
type of stress a component is subject to, different models are 
used to simulate crack growth. For stresses that are approxi- 
mately constant in time, creep crack growth models are applied. 

A typical creep crack growth model is given by — - = y(C,y", where 

dt 

a is the crack length, Ct a so-called crack tip parameter that 
depends on the component geometry and applied stress, y ^ mate- 
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rial creep constant, and m ~ n/(n+l) another constant with n 
being the Norton exponential constant . For stresses that are at 
least approximately cyclic, fatigue crack growth models are 
used. A typically used fatigue crack growth model is given by 

5 — = — — , where AX is the amplitude of the 

^crit 1 
K 

max 

stress cycle, a the crack length and N the number of cycles and 
the remaining variables are component specific constants. 

The respective models are well-known to a person skilled in the 
art of mechanical engineering, and are e.g. described in the 

10 books "^Fracture Mechanics: Fiindamentals and Applications", by 
T. Anderson, CRC Press, 2nd ed. , Cleveland, Ohio, 1995; '*Ad- 
vanced Fracture Mechanics" by M. Kanninen and C. Popelar, Ox- 
ford University Press, New York, 1985; ^Non- linear Fracture Me- 
chanics for Engineers", by A. Saxena, CRC Press, Cleveland, 

15 Ohio, 1998; and ^Damage Mechanics and Life Assessment of High- 
Temperature Components" by R. Vishwanathan, ASM International, 
Ohio, 1989. Implementation of the modelling methodology into a 
crack-growth simulator is again straightforward to a person 
skilled in the art of computer programming. 

20 To further illustrate how the crack growth simulations may be 
carried out. Fig. 1 shows a potential location and geometry of 
a crack 3 that has been growing for some time in a turbine 
blade 1, said turbine blade 1 being fixed to a rotor 2. When 
the turbine is in operation, the blade 1 is subject to a ten- 

25 sion P and a bending moment if, which are indicated by arrows. 
Fig. 2 shows a cross section of a pipe 4, with an example of a 
crack growth geometry for a crack 5 that has been growing for 
some time in a radial direction of the pipe 4. The pipe 4 is 
sxibject to an internal pressure p exerted by flow of gas 

30 through the pipe 4 . 
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The life expectancy of the component obtained from the crack- 
growth simulation is then compared to a given, desired value. 
If the life expectancy of the component obtained from the 
crack-growth simulation is shorter than the desired value, the 
structure of the component is modified, and the computation of 
stress as described in the second step above and the determina- 
tion of the life expectancy by simulation of crack-growth as 
described in the third step above, respectively, are repeated 
with the new, modified structure of the component. This process 
is continued in an iterative manner until the life expectancy 
exceeds the desired value. 

If the life expectancy is longer than the desired value, the 
component can be used as is, i.e. the structure of the compo- 
nent does not require further modification. This requirement is 
already fulfilled after the 0-th iteration if the life expec- 
tancy obtained with the initial structure of the component is 
longer than the desired value. In such cases, the initial 
structure need not be modified. 

In a preferred variation of the method, the structure is param- 
eterized by means of at least one structural parameter, pref- 
erably a plurality of structural parameters. Structural parame- 
ters may refer to an inner or outer diameter of a tube, a wall 
thickness of a casing, a blade/vane height, width, or cord, a 
rotor mass, structural interconnections, a material composi- 
tion, in particular a ratio of a certain sxibstance in a compos- 
ite or an alloy, a material tempering, etc. In the iterative 
process, the structure of the component is preferably modified 
by modifying one or more of the structural parameters. 

In another preferred variation of the method, instead of com- 
paring the life expectancy to a given, desired value, it is 
checked whether the life expectancy falls into a given inter- 
val. By applying this variation of the method to a number of 
components that constitute a subunit of the steam turbine, it 
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can be assured that all components in this subxinit have the 
same life expectancy. This is of particular advantage for sub- 
xmits that are regularly replaced, as it ensures that no compo- 
nents with a large remaining life expectancy are replaced un- 
5 duly if all components in the subunit are provided with 
approximately the same life expectancy, i.e. within an interval 
that is small compared to the life expectancy. 

In another preferred variation of the method, not only the com- 
putation of stress and the determination of the life expectancy 

xo are repeated after each modification of the structure of the 
component in each iteration, but also the process variables are 
re -computed by means of the process simulator whenever the 
structure of the component has changed significantly since a 
previous computation of the process variables. If the structure 

15 is parameterized, a significant change can advantageously be 
defined by specifying for each parameter an allowable relative 
change, beyond which a corresponding change in the structure 
has to be regarded as significant. 

In another preferred variation of the method, the process vari- 
20 ables are re-computed by means of the process simulator in each 
step of the iteration, i.e. after each modification of the 
structure of the component . 

A computer program product according to the invention comprises 
a computer readable medium, having thereon computer program 
25 code means that, when loaded onto a computer, make said com- 
puter execute the method or any of the preferred variations as 
described above in this section. 
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List o£ Reference Signs 

1 Turbine blade 

2 Turbine rotor 

3 Crack in turbine blade 
P Tension 

M Bending moment 

4 Pipe 

5 Crack in pipe 

p Internal pressure . 
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